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Abstract. We discuss an idea for how accreting millisecond pulsars could contribute to the under- 
standing of the QCD phase transition in the high-density nuclear matter equation of state (EoS). 
It is based on two ingredients, the first one being a "phase diagram" of rapidly rotating compact 
star configurations in the plane of spin frequency and mass, determined with state-of-the-art hybrid 
equations of state, allowing for a transition to color superconducting quark matter. The second is the 
study of spin-up and accretion evolution in this phase diagram. We show that the quark matter phase 
transition leads to a characteristic line in the £2 — M plane, the phase border between neutron stars 
and hybrid stars with a quark matter core. Along this line a change in the pulsar's moment of inertia 
entails a waiting point phenomenon in the accreting millisecond X-ray pulsar (AMXP) evolution: 
most of these objects should therefore be found along the phase border in the Q. — M plane, which 
may be viewed as the AMXP analog of the main sequence in the Hertzsprung-Russell diagram for 
normal stars. In order to prove the existence of a high-density phase transition in the cores of com- 
pact stars we need population statistics for AMXP's with sufficiently accurate determination of their 
masses and spin frequencies. 
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INTRODUCTION 

Accreting compact stars in low-mass binary systems undergo a stage with disc accretion 
leading to both spin-up and mass increase. Initial indications for spin frequency clus- 
tering in low-mass X-ray binary (LMXB) systems, reported by measurements with the 
Rossi-XTE, have lead to the interpretation of such a correlation as a waiting-point phe- 
nomenon where stellar configurations cross the border between pure neutron stars and 
hybrid stars in the spin frequency-mass plane A systematic analysis of the critical 
line for a deconfinement phase transition in the "phase diagram" for accreting compact 
stars [0] has revealed that the suggested population clustering due to the phase transi- 
tion will instead lead to mass clustering. For homogeneous interiors like in the case of 
strange stars, however, such an effect shall be absent id]. For generic polytropic forms 
of the EoS of quark and hadronic matter, the relationship between softness or hardness 
of the EoS and the structure of this phase diagram has been demonstrated in Ref. [4]. 

In the present contribution we give an updated view on these ideas, based on a recently 
developed hybrid EoS JH] which fulfills constraints from compact star observations and 



heavy-ion collision (HIC) experiments This example demonstrates the compatibility 
of an onset of deconfinement in stars with typical masses ~ 1.4 Mq with a high 
maximum mass ~ 2.0 Mq. We speculate that the class of objects for which mass 
clustering due to a phase transition in their interior applies might include also most 
of the first-bom neutron stars in low-kick, non-excentric double neutron star as well as 
in pulsar- white dwarf systems which have undergone a mass accretion stage f/l]. 

Upon further elaboration and proper selection of the pulsar population, the mass 
clustering phenomenon could prove to be a direct observation of a phase transition in 
compressed nuclear matter, such as the QCD chiral symmetry restoration transition. 



HYBRID EOS: MASS-RADIUS VS. FLOW CONSTRAINTS 

One of the most challenging tasks in fundamental nuclear and particle physics is 
the delineation of the border between hadronic matter and quark-gluon matter in the 
temperature-density plane: the QCD phase diagram. At high temperatures and low 
baryon densities, numerical simulations of QCD as a lattice gauge theory indicate 
that both, chiral symmetry restoration and deconfinement are crossover transitions and 
their critical temperatures coincide at a value of T^, — T^j = 196 MeV, obtained by the 
Bielefeld-Brookhaven-Columbia-Riken collaboration This value is consistent with 
a statistical model analysis of the freeze-out temperature 7/ = 160 MeV from hadron 
production in nuclear collisions at RHIC for the highest presently available c.m.s. en- 
ergy of = 200 Ge V M, [H • 



At finite densities and T ^ 0, the situation is quite different. Both lattice QCD 
simulations and heavy-ion collision experiments cannot access this region and many 
questions are yet unanswered, like: 

(i) Is (are) the transition(s) of first order so that there must be a critical endpoint in the 



QCD phase diagram [|12|1? 



(ii) Is quark matter at low T a color superconductor and how does such a property 



manifest itself 111311 ? 
Other points of discussion which we will not discuss in detail here are: 
(iii) Do all quark flavors appear simultaneously or rather sequentially [ 14l [l5l [l6Lfl71? 



(iv) Do chiral and deconfinement transition happen at the same critical chemical poten 



tial jl), = jld or is there a so-called quarkyonic phase [11811 ? 

At present, the description of cold, dense matter phases within QCD is out of reach. 
Therefore, effective models for quark matter like those of the Nambu-Jona-Lasinio 
(NIL) type [19] may help to obtain quantitative estimates for the dense matter EoS and 
phase transitions, once the free parameters can be fixed. In this situation, astrophysical 
observations of compact stars may provide constraints for the behavior of matter under 



We want to remark that a considerably lower value of = 151 MeV has been reported by Fodor et al. 
i^. The discrepancy to Cheng et al.|@] is not yet resolved. 



high compression and isospin asymmetry, complementary to data from lattice QCD and 
heavy-ion collisions. 
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FIGURE 1. Left: Hybrid EoS with DBHF nuclear matter and color superconducting NJL quark matter 
fitted to obey the flow constraint |22]. Right: Mass-Radius constraints from thermal radiation of the 
isolated neutron star RX J1856.5-3754 (grey hatched region) and from QPOs in the LMXB's 4U 0614+09 
(green hatched area, the "wedge") and 4U 1636-536 (orange hatched region) which shall be regarded as 
separate conditions to the EoS, see Ref. [5]. The controversial interpretation that the 1122 Hz pulsations 
for XTE J1739-285 correspond to its spin frequency |23] entails a stringent mass-radius constraint for 
this compact star |24i] . The theoretical M — R relations correspond to hybrid star configurations with a 
DBHF hadronic shell and a color superconducting NJL quark matter core. The underlying hybrid EoS 
is the isospin-asymmetric generalization of that in the left panel, under )3 -equilibrium with electrons and 
muons. Shown are also lines of constant surface redshift z- 

In the spirit of two-phase models for the dense matter EoS, the nuclear matter EoS can 
be described by the ab-initio Dirac-Brueckner-Hartree-Fock (DBHF) approach [|20| . I21I1 
using the Bonn-A nucleon-nucleon potential and the quark matter EoS is given by a 
color superconducting three-flavor NJL model [15] augmented by a vector meanfield 
contribution [5] which stiffens the deconfined phase. In the left panel of Fig. [U we illus- 
trate how the coupling strengths in the scalar diquark (rjo) and in the vector meson {r]v) 
channels as free parameters of the NJL model description can be fixed by the require- 
ment that the isospin-symmetric hybrid EoS resulting from a Maxwell construction with 
the DBHF EoS is as stiff as possible but still compatible with the constraint region from 
HIC flow data at AGS and SIS energies [22]. Resulting preferable parameter values are 
rjY = 0.5 and rjo in the range rjD = 1.02... 1.03. The onset of the hadron-to-quark mat- 
ter transition in symmetric nuclear matter is at a density of = 0.55 fm^^ = 3.5 hq, 
where no = 0. 16 fm^-^ denotes the nuclear matter saturation density. The generalization 
of this hybrid EoS to the isospin-asymmetric case under /3 -equilibrium with electrons 
and muons can be used to predict the critical mass of a neutron star for the onset of 
deconfinement in its interior in the range of Merit ~ 1.35 ... 1.0 M©, see also [5]. The 
corresponding critical density for the hadron-to-quark matter transition in compact stars 
is at about Ucnt = 0.4 fm^-^ = 2.5 uq. In the right panel of Fig. [H we show the sequences 
of hadronic compact stars (red line) and the hybrid star sequences with color supercon- 
ducting quark matter cores (black and blue lines). We would like to stress two results: 
(i) we predict that compact stars in the typical mass range of 1.35 ±0.1 Mq may have 
quark matter interiors, and (ii) the maximum mass for compact stars with quark matter 



interiors exceeds 2 Mq. The latter result disproves the claim (l^ that the occurrence of 
deconfined quark matter in compact stars would necessarily lead to a softening of the 
EoS which excludes Jhe possibility of hybrid star configurations with large masses and 
radii, see alsoRef. i26ll. 

However, this demonstrates only that the possibility of compact stars having quark 
matter interiors cannot be excluded by the measurement of a large mass. What does 
it take to discover a high-density phase transition by astrophysical observations? With 
the newly developed hybrid star EoS at hand we want to revisit the suggestion that 
a clustering of frequencies [T] and/or masses in the population of compact stars 
in LMXB's may be seen as a signal for deconfinement and add a new twist to this 
hypothesis by suggesting to include double neutron stars into these considerations! 

POPULATION CLUSTERING AS A SIGNAL FOR 
DECONFINEMENT 

Our investigation is based on the classification of rotating compact star configurations 
in the plane of angular velocity Q. and baryon number (or gravitational mass M), the 
so called phase diagram for compact stars [27]. It is defined by at least three lines: the 
maximum frequency ^max{N) for which stable rotation without mass shedding can be 
sustained, the maximum baryon number A^max(^^) the star can carry without undergoing 
gravitational collapse and the critical line Ncnt{^), which separates the region of quark 
core configurations from the hadronic ones, see Fig. [2l It has been shown that the latter 
line is correlated with the local maxima of the moment of inertia with respect to changes 
of the baryon number at given Q. due to the change of the internal structure of the 
compact object at the deconfinement phase transition. Therefore, we expect that the 
rotational behavior of these objects changes in a characteristic way when this line is 
crossed. Early suggestions of a deconfinement signal following from this characteristic 
behavior have considered the spin-down of isolated radio pulsars without mass accretion 
for which a characteristic deviation of the braking index from the value n = 3 for dipole 
emission shall signal the transition [28, 29]. It has also been noted that the star has to 
spend about 10^ yr for crossing the configuration border in the phase diagram, the typical 
time it takes the star to loose by dipole emission the amount of angular momentum 
A7 = QA/ which corresponds to the change in the moment of inertia AI sue to the phase 
transition in the star's interior. The consequence will be an increase of the population 
of stars at this critical line which could be observed. Provided that a sufficiently large 
number of accretors will be discovered and their masses and spin frequencies [2] could 
be determined. Then the phase transition would reveal itself by a population clustering 
along the line A^crit(^) in the phase diagram, which according to the above results for 
nonrotating stars should be well separated from the black hole limit Njnaxi^)- 

The main problem with this signal is the shape of the critical line in the phase diagram 
which for typical hybrid EoS [4] would suggest a mass clustering [2] rather than a 
frequency clustering [1]. In order to measure the star mass, however, one would need 
a companion star which for isolated pulsars is absent! 

Therefore, the suitable population for which this statistical phase transition test shall 
be applicable are compact stars in binary systems with mass transfer via Roche-lobe 



overflow. There are two such systems which we want to focus on: (i) LMXB's [|30ll and 
(ii) double neutron stars (DNS's), see [7] for a recent discussion. 

The spin evolution of a compact star under mass accretion from a low-mass compan- 
ion star can be regarded as a sequence of stationary states of configurations (points) in 
the phase diagram. It is governed by the change in angular momentum of the star 

j^{I{N,a)a)=K,^,, K,u=VGMM^ro-Noni, (1) 

where ^ext denotes the external torque due to both the specific angular momentum 
transfered by the accreting plasma and the magnetic plus viscous stress given by A^out = 
KjA^r^^, fc = 1/3 [31]. For a star with radius R and magnetic field strength B, the 

1 /3 

magnetic moment is given hy il = B and = [GM/Q?) is the co-rotating radius, 
see iliB] and references therein for details. From Eq. ([T]) follows the evolution equation 
for the angular velocity 



dt ~ IiN,Q.)+a{dIiN,a)/dQ.)r, ' '^'^'^'^^''> ^ J^' 

Solutions of (l2l) are trajectories in the Q — A'^ plane describing the spin evolution of ac- 
creting compact stars. Since I{N, Q.) exhibits characteristic functional dependences fr^ 
at the deconfinement phase transition line A^crit(^) we expect observable consequences 
in the P — P plane when this line is crossed. 

In our model calculations we assume that both the mass accretion and the angular 
momentum transfer processes are slow enough to justify the assumption of quasistation- 
ary rigid rotation without convection. The moment of inertia of the rotating star can be 
defined as I{N, Q.) = J{N, , where J{N, ^) is the angular momentum of the star. 



For a more detailed description of the method and analytic results we refer to [1290 and 
references therein. The time dependence of the baryon number for the constant accreting 
rate A^^ is given by N{t) = N{to) + {t — to)N and for the magnetic field of the accretors 
we consider the exponential decay B{t) = [5(0) — 5oo] exp(— r/Tg) + B^. We solve the 
equation for the spin-up evolution ^ of the accreting star for decay times Ts = 10^ yr 
and initial magnetic fields in the range 0.2 < 5(0) [TG] < 4.0. The remnant magnetic 
field is chosen to be 5., = 10-4tg[32], where 1 TG= lO^^ G. 

The question arises whether there is a characteristic feature in the spin evolution 
when the trajectories traverse the critical phase transition line. In order to perform a 
more quantitative discussion of possible signals of the deconfinement phase transition 
we investigate the spin-up evolution for stars with A^(0) = 1.4 Nq and Q.{0) = 1 Hz in 
the initial state. In the case of high accretion rate {N = 10^^A/^3/yr, e.g. for Z sources) 
and long-lived magnetic field (Tg = 10^ yr) there is a dip in the period derivative when 
the star evolves into the quark core region of the phase diagram. This feature can be 
quantified by the distribution of a waiting time % = \P/P\ =^/CL 'm the Q — A'^ plane. 

In Fig. |2]we show contours of waiting time regions in the phase diagram. The region 
of longest waiting times is located in a narrow branch around the phase transition border 
and does not depend on the evolution scenario after the passage of the border, when the 
depopulation occurs and the probability to find an accreting compact star is reduced. 
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FIGURE 2. Regions of waiting times in the phase diagram for compact hybrid stars for a magnetic 
field decay time of = 10^ yr and a constant mass accretion rate of = 10^** A^o/yr- For an estimate 
of a population statistics we show the region of evolutionary tracks when the interval of initial magnetic 
field values is restricted to 0.6 < B(0)[TG] < 1.0. Note that the probability of finding a compact star in 
the phase diagram is enhanced in the vicinity of the critical line for the deconfinement phase transition 
^crit(^i) by at least a factor of two relative to all other regions in the phase diagram. The correlation of the 
distribution of objects with the line for a deconfinement phase transition could be interpreted as a waiting 
point phenomenon during the accretion evolution of the compact stars (mass clustering), see L2J. 



Another smaller increase of the waiting time and thus a population clustering could 
occur in a region where the accretor is already a quark core star. For an estimate of the 
magnetic field influence we show in Fig. [2] also the region of evolutionary tracks when 
the values of initial magnetic field vary within 0.6 < 5(0) [TG] < 1.0. 

As a strategy of search for QCSs we suggest to select from the LMXBs exhibiting the 
QPO phenomenon those accreting close to the Eddington limit [|33l] and to determine 
simultaneously the spin frequency and the mass [|34l] for sufficiently many of these 
objects. The emerging statistics of accreting compact stars should then exhibit the 
population clustering shown in Fig. [21 when a deconfinement transition is possible. If 
a structureless distribution of objects in the Q. — N plane will be observed, then no 
firm conclusion about quark core formation in compact stars can be made as, e.g., for 
strange stars [3]. There is a problem with the simultaneous analysis of masses and spin 
frequencies for LMXB's which is still rather dependent on the model employed [35, 36]. 

Following van den Heuvel Ivj], we may suggest to apply our approach also to the 
population of the first-born stars in double neutron star systems with low excentricity 
and pulsar- white dwarf double systems, as they should have undergone a mass accretion 
stage similar to the one described above which is a prerequisite for the applicability 
of the population clustering signal of deconfinement. For the specific example of the 
hybrid EoS discussed above with a DBHF hadronic phase and a color superconducting 
stiff quark matter phase, we show in Fig. [3] the phase diagram of rotating compact stars 
together with the masses and spin frequencies obtained for double neutron star systems. 
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FIGURE 3. "Phase diagram" for compact stars in the plane of rotation frequency D. vs. mass M. 
Stable stars are found in the region bordered by the maximum rotation frequency Q.max{M) for which 
mass-shedding from the star equator occurs and the maximum mass Mniax(ii) for which the gravitational 
instability against collapse to a black hole sets in. The dash-dotted (dash-double-dotted) line marks the 
onset of a quark matter core in the stars interior for a diquark coupling parameter rjo — 1 .017 (rjo — 1 -02), 
see[T]for the corresponding mass-radius relations of the nonrotating star sequences. The data with error 
bars con^espond to binary radio pulsars (black dots) [7] and neutron stars with a white dwarf binary 
(green squares) ll24ll . The correlation of the distribution of objects with the line for a deconfinement 
phase transition could be interpreted as a waiting point phenomenon during the accretion evolution of the 
compact stars (mass clustering), see 



We observe an interesting correlation of the distribution of these objects with the critical 
deconfinement phase transition line for diquark coupling in the range r}o= 1.017...1.02. 
There might be other reasons for the mass clustering of spinning pulsars [7], but the 
suggestion to relate it to a phase transition in the compact star interior can not be 
excluded. 



SUMMARY 

The model independent result of our study is that a population clustering in the phase 
diagram for accreting compact stars shall measure the critical line A'crit(^^) which sep- 
arates neutron stars from hybrid stars where the shape of this curve can discriminate 
between different models of the nuclear EoS at high densities. 

For the new hybrid equation of state discussed in this contribution, we expect the 
suggested population clustering as a signal of the deconfinement transition in the mass 
range M ~ 1 .3 ... 1 .4 Mq which might therefore serve as one aspect for the explanation 
of the well-known mass clustering in double neutron stars in that same mass range. The 
above rather fresh ideas suggest that future observational programs may contribute to 
unraveling most actual problems of fundamental physics of dense baryonic matter. 
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